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Abstract. We study the quantum analog of stirring of water inside a cup using a spoon. This can be regarded 
as a prototype example for quantum pumping in closed devices. The current in the device is induced by 
translating a scatterer. Its calculation is done using the Kubo formula approach. The transported charge is 
expressed as a line integral that encircles chains of Dirac monopoles. For simple systems the results turn 
out to be counter intuitive: e.g. as we move a small scatterer "forward" the current is induced "backwards". 
One should realize that the route towards quantum-classical correspondence has to do with "quantum chaos" 
considerations, and hence assumes greater complexity of the device. We also point out the relation to the 
familiar S matrix formalism which is used to analyze quantum pumping in open geometries. 



1. Introduction 

Consider a closed ring that contains particles (Fig. la). Assume that one wants to create a 
current in this ring. If the particles are charged then one way to do it is by creating an electro 
motive force (EMF). This can be induced by varying an Aharonov-Bohm flux <t>, such that by 
Faraday's law EMF = -<t>. But there is another way to create a current that does not involve 
EMF, and hence does not assume charged particles. The idea is to change in time the scalar 
potential V(r;Xi(f),X2(tJ). Here r is the coordinate of a representative particle in the ring, 
while X\ and X2 are some control parameters. By making a cycle in the (Xi , X2) space we can 
push non-zero net charge Q through the system. Thus an "AC driving" gives rise to a "DC" 
component in the current. This is known in the literature as "quantum pumping". 




Fig.l. Models for the analysis of quantum stirring, (a) Upper panel: A scatterer (big black dot) is translated inside a 
Sinai billiard. A chaotic trajectory of a representative particle in this billiard is illustrated, (b) Lower panels: Network 
models for quantum stirring. The scatterer (big black dot) is translated along one of the bonds. The vertical dotted 
line is the section through which the current is measured. From left to right: chaotic network; double barrier model; 
triple barrier model. 
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Fig.2. (a) Left panel: A schematic representation of the network model. The vertical dotted line is the section through 
which the current is measured. The moving scatterer is indicated by its transmission go, while X[ is its displacement 
along the bond, (b) Right panel: The corresponding open geometry where the left and right leads are connected to 
reservoirs with the same chemical potentials. 

In this paper we would like to consider a prototype pumping problem, which we call 
"quantum stirring". It is the simplest scheme to create a current with a non-vanishing DC 
component. Referring to Fig.2 we define X\ as the location of a scatterer, while X2 is its 
"size". By "size" we mean either the cross section or the reflection coefficient. One can 
regard the scatterer as a "piston" or as a "spoon" with which it is possible to "push" the 
particles. A prototype example for a pumping cycle is illustrated in Fig.3. During the main 
stage of the cycle the scatterer is translated to the right a distance AX\ . Consequently a charge 
Q is transported. In the second stage the size of the scatterer is "lowered", and it is displaced 
back to its original location, where its original "size" is restored. By repeating this cycle many 
times we can create a current with a DC component. 

In the following analysis we assume that the system consists of non-interacting spinless 
particles. All the particles have (formally) charge e, even if they are not actually charged. We 
assume that there is no magnetic field in the system. Still, for the sake of a later mathematical 
formulation, it is convenient to introduce a third parameter X? = O, where cD is an Aharonv- 
Bohm flux. The pumping cycle in the (X\, X2, X3) space is illustrated in Fig.3. 




Fig.3. A prototype example for a pumping cycle. During the main stage of the cycle the scatterer is translated to the 
right a distance AXi . Consequently a charge Q is transported, (a) Left panel: The pumping cycle in the 2-dimensional 
(Xi , X2) plane, (b) Right panel: The same pumping cycle in the three dimensional (Xi , Xj, X3) space, where X3 = <D 
is the Aharonov Bohm flux via the ring. 
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Fig.4. (a) Left panel: The calculation of the charge Q is a line integral over G that can be regarded as a calculation 
of the flux of B via a two dimensional curve, ds is a normal vector to the pumping cycle. The black dot in 
the middle symbolizes the presence of "magnetic charge" which is characterized by a density (t(Xi,X2). In the 
quantum mechanical analysis this should be understood as the density of "Dirac chains", (b) Right panel: In the 
embedding (X\,X?,X)) space the magnetic charge is organized as vertical charged chains. Each chain consists of 
"Dirac monopoles" which are located at X" points where an occupied level has a degeneracy with a nearby level. The 
ellipse represents a possible pumping cycle that may encircle either one or many chains. 



1.1. Linear response theory and the Dirac chains picture 

We are going to analyze the stirring problem within the framework of linear response theory. 
If we have EMF then we expect to get in the DC limit Ohm law I = -G4>, while if we change 
slowly either X\ or X 2 we expect to get in the DC limit I - -G 1 X\ or I — -G 2 X 2 respectively. 
So in general we can write 

Q = (j) Idt=- (^(G^dXy + G 2 dX 2 ) = (j)B -£s = jj o-(X l ,X 2 )dX l dX 2 (1) 

In the second expression we define the normal vector ds = (dX 2 , —dX\) and use the notations 
B\ = -G 2 and B 2 — G 1 . See Fig. 4a for an illustration. The third expression is obtained via 
the two dimensional version of the divergence theorem. If we regard B as a fictitious magnetic 
field, then cr is the two dimensional density of magnetic charge. 

It turns out that in the strict adiabatic limit the vector field B is related to the theory of 
Berry phase (DEI- The formulation of this relation is as follows. Assume that the system is 
adiabatically cycled in the (X\ , X2, X3) space. In such case the Berry phase can be calculated as 
a line integral over a "vector potential" (also called "1 form") A. This can be converted by the 
Stokes theorem into a surface integral over a "magnetic field" (also called "2 form") B. The 
B field is defined as the "rotor" of A. It is a divergence-less field but it can have singularities 
which are known as "Dirac monopoles". These monopoles are located at X points where an 
occupied energy level has a degeneracy with a nearby level. Because of O i-> O + (2nh/e) 
gauge invariance the Dirac monopoles form vertical chains as illustrated in Fig. 4b. Hence we 
have a distribution of what we call "Dirac chains" 00), which is characterized by a density 
cr(X u X 2 ). 
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1.2. Background and objectives 

Most of the literature about quantum pumping deals with the open geometry of Fig. 2b. The 
most popular approach is the S -scattering formalism which leads to the Biittiker, Pretre and 
Thomas (BPT) formula |5 6| for the generalized conductance G. The BPT formula, is 
essentially a generalization of the Landauer formula. In previous publications 1 7 4 1 we have 
demonstrated that the BPT formula can be regarded as a special limit of the Kubo formula. 
Our Kubo formula approach to pumping 1 3 8 1 leads to "level by level" understanding of the 
pumping process, and allows to incorporate easily non-adiabatic and environmental effects. 
In the strict adiabatic limit it reduces in a transparent way to the theory of adiabatic transport 
j9]IT0l. also known as "geometric magnetism |2|. On the other hand, in the non-adiabatic(I) 
"DC limit" of an open geometry it reduces to the S matrix picture, hence resolving some 
puzzles that had emerged in older publications. 

The question "how much charge is pushed by translating a scatterer" has been addressed 
in Ref. 1 1 1 1 in the case of an open geometry using the BPT formula. We have addressed the 
corresponding problem of quantum stirring in closed geometry in a previous short publication 
1121 . but the connection with the Dirac chains picture has not been illuminated. Furthermore, 
in | only the quantum chaos limit was considered. 

In the present publication we put an emphasis on clarifying the route towards 
quantum-classical correspondence (QCC). We shall see that quantum mechanical effects are 
pronounced in simple systems. As the system becomes more chaotic QCC emerges. The 
Dirac chains picture leads to new insights regarding the route towards QCC. These insights 
are easily missed if we stick to the formal Green function calculation of our earlier work 1121 . 
From the above it should be clear that the main objectives of the present study are: 

• Derivation of a classical formula for Q (assuming a stochastic picture). 

• Derivation of a quantum result for Q using the Dirac chains picture. 

• Exposing some counter-intuitive results for Q in the case of the simplest models. 

• Illuminating the route towards QCC as we go from "simple" to "chaotic" systems. 
We note that in 1 12 1 we have presented the classical formula for Q without the derivation. 

1.3. Physical motivation and experimental feasibility 

In the previous section we have explained the theoretical motivations for dealing with the 
stirring problem. In the present section we would like to further discuss the practicality of this 
line of study, and the feasibility of actual experiments. 

It is quite clear that the main focus of today's experiments is on open devices (with 
leads), whereas our interest is in closed devices. Our believe is that "wireless" mesoscopic 
or molecular size devices are going to be important building blocks of future "quantum 
electronics". This is of course a vision that people may doubt. However, on the scientific 
side our task is to analyze its feasibility. 

It is possible to fabricate closed mesoscopic rings, and to measure the persistent or the 
induced currents. Experiments with closed devices have been performed already 10 years 
ago. As an example we mention Ref.| 13 1 where a large array of rings has been fabricated. 
The current measurement has been achieved by coupling the rings to a highly sensitive 
electromagnetic superconducting micro-resonator. 

The conceptually simplest way to drive a current is by inducing an electro motive force 
(EMF). In the setup of Ref.| 13 1 the EMF has been induced by a "wire" that spirals on top of 
the array. In our view an attractive alternative option would be to induce currents by changing 
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gate voltages so as to induce stirring. The advantage of such a possibility for the purpose of 
integrating wireless devices in future quantum electronics is quite obvious: It is much easier 
to control gate voltages than fluxes of magnetic field. 

As far as electronic devices are concerned there is no question about the feasibility 
of realizing quantum stirring by manipulating gate voltages, and measuring the electrical 
currents. But we would like to argue that such possibility is open also in case of neutral atoms. 
It is well known that "billiards" that confine cold atoms can be realized and manipulated 
1 141 1 1 51 . Furthermore, there is no question regarding the possibility of creating a "moving" 
optical barrier so as to create a stirring effect. There are variety of techniques to measure the 
induced neutral currents. For example one can exploit the Doppler effect at the perpendicular 
direction, which is known as the rotational frequency shift 1 16 1. 

There is one more issue which might be of relevance in case of an actual experiment. The 
Kubo formalism assumes that the system settles into a steady state, whereas the preparation 
in case of an actual experiment is not very well controlled. We would like to argue that the 
results of the linear response analysis are quite robust. This issue is discussed in section 4 of 
Ref.| 17 1 : What we get for Q in the Kubo analysis is not merely a formal result, but rather a 
prediction that has an actual physical significance. 

1.4. Outline 

In the first part of this paper we review the result for G in the case of an open system using 
the BPT formula. Then we present two equivalent derivations of the corresponding classical 
result in the case of a closed geometry. We use the term "classical" in the Boltzmann sense. 
This means that interference within the ring is neglected, while the reflection by the scatterers 
("cross section") is calculated quantum mechanically. The first derivation is based on a direct 
solution of a master equation, while the second is a straightforward application of the Kubo 
formula. The classical calculation implies an expression for the density o~(X\,X2) of the 
monopoles. The BPT formula implies cr{X\,X2) that can be regarded as a special case of this 
calculation. 

In the second part of this paper we turn to the quantum mechanical analysis. As a 
preliminary stage we discuss the general conditions for having a degeneracy point X in the 
case of a one dimensional ring. Then we review how the pumped charge Q can be estimated 
by calculating a line integral that encircles "Dirac chains". Thus we realize that we have to 
figure out what' the distribution o-{X\ , Xi) of these chains looks like. Specifically, we consider 
the model systems that are illustrated in Fig.l, and schematically in Fig. 2. The simplest is a 
ring where both gi and go are modeled as delta barriers. The result for Q is quite remote from 
the classical expectation. Consequently we try to figure out what happens to cr{X\ , X2) as the 
system becomes more complex: First we add a second fixed barrier, and finally we consider 
what happens in the case of a "chaotic" barrier which is modeled using random matrix theory. 
We make it clear that the route to the classical limit is intimately related to so called "quantum 
chaos" considerations. 
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2. Pushing particles in an open geometry 

Let us consider the model of Fig. 2b, where we have a scatterer within a single mode wire 
which is connected to two reservoirs with the same chemical potential. In this section 
we assume non-interacting spinless electrons and zero temperature Fermi occupation. The 
scatterer is described by 

V(r,X 1 ,X 2 )=X 2 8(r-X{) (2) 

Hence, for some fixed values of X\ and X2 its transmission is 

(3) 



8o(X 2 ) = 



1+1 



where m is the mass of the particle and k F is the Fermi momentum. From now on we work 
with units such that H = 1 . The S matrix of the scattering region can be written in the general 
form 



•■'1 V <#jH (4) 

where y is the total phase shift, a is the reflection phase shift, and <p = e<S)/H represents the 
flux which we assume to be zero. In the setup of Fig. 2b the length of the right lead is La - X\ 
and the length of the left lead is Lb +X\. Hence 

8 = 80 (5) 

y = k F (L A +L f f)-arctan(-^-X 2 ) (6) 



\H 2 k F 

a = k(L A - L B ) -UcpXi (7) 

Now that we know the dependence of the S matrix on the parameters (Xi , X 2 ), the calculation 
of G is quite straightforward. We use the BPT formula 

G;=-^frace(p lead |V] (8) 



2ni \ ™dXj 

where projects on the channels of the lead where the current is measured. As indicated in 
Fig. 2b the current is measured via a section which is located on the right lead. Using the BPT 
formula we get 

G l = -(l-8o)-h (9) 

° 2 = -*°4^ (10) 

where v F is the Fermi velocity corresponding to k F . The result for G 1 is our main interest. 
It has been discussed in Ref. Bl II . where the term "snow plow" has been coined in order 
to describe its physical interpretation. Namely, for zero temperature Fermi occupation the 
density of electrons in the wire is k F /n. Therefore the number of electrons that are pushed 
by the scatterer is dN = (k F /n) X dX\ . If the transmission of the scatterer is not zero, some 
of the electrons pass through it and consequently we have to multiply dN by the reflection 
probability 1 - go. 
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3. Stirring of particles in a closed geometry 



Let us consider the model of Fig. 2a, where the system is closed. We assume that the 
transmission of the ring without the moving scatterer is gf, while the transmission of the 
scatterer itself is go. In the following two subsections we shall present two optional derivations 
of the "classical" result for G. We use the term "classical" in the Boltzmann sense. Namely, 
we regard the scattering from either gf or go as a stochastic process. Thus interference within 
the arms of the ring is not taken into account. For sake of comparison with the BPT-based 
result we still assume zero temperature Fermi occupation (while in later sections we shall 
allow any arbitrary occupation). Within this framework we obtain: 



G - - 



G - - 



(1 



go)gf 



go + g\ 
(1- 



d -2g gf 



n 



?[ l )go 



go + gf - 2g g c , 



Anhv ¥ 



(11) 



(12) 



We note that the amount of charge which is pushed by translating a scatterer a distance AX\ 
can also be written as 1121 



Q 



-G l AX 



1 



-go 



So 



gT 



1 



■gT 



71 



(13) 



1 ~gT 




1 - go 




I 1 -*?] 


gT 




go 


+ 


[ gf \ 



where gj is the overall transmission of the ring (including the moving scatterer) if it were 
opened: 



(14) 



As expected the charge Q which is transported as a result of an X\ displacement depends in a 
monotonic way on the reflection coefficient 1 - go- It monotonically increases from zero, and 
attains half of its maximal value for go = gf. A plot of Q versus the "size" of the scatterer is 
presented in Fig. 5 for three representative values of gf. We also plot Q against X2, assuming 
that the scatterer is modeled as a delta function. 




Fig.5. Plots of Q as a function of the "size" of the scatterer. We use arbitrary units such that Q = 1 in the maximum, 
(a) Left panel: Q is plotted against the reflection coefficient (1 - go) for gf = 0.1, for gf = 0.5, and for gf = 0.9. 
The dotted lines highlight that Q for go = gf is half its maximum value. Note that the BPT based result corresponds 
to gf = 0.5. (b) Right panel: Here Q is plotted against Xi assuming that the scatterer is a delta function, and 
setting rr\/(h 2 k F ) = 1. 
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It is important to realize that the result for an open geometry is formally a special case 
corresponding to gf = 1/2. This value of gf means that memory is completely lost once a 
particle is scattered by the "surroundings". Namely, if gf = 1/2 then after a collision a particle 
has equal probability to go in either direction, and any information about its initial direction 
is lost. This observation generalizes our discussion in Ref.| 18 1 regarding the relation between 
the Kubo and the Landauer conductance. 

The classical expression for G implies the following result for the density cr{X\,X2), 
which is illustrated in Fig. 6. 

2(1-*?)*? 



<r(X 1 ,X 2 ) = 



dB 2 



em 



dX 2 nh 1 



1 + 



((ft*) 2 ->); 



m 



■X, 



(15) 




Fig.6. The classically deduced density cr as a function of X2 for gj =0.1, for gj = 0.5, and for gj = 0.9. We use 
arbitrary units for o", and set m/(H 2 k F ) = 1. The dotted vertical lines correspond to the median X2 values which are 
determined by the equation go(X2) = gf. 

In the following sections we give two optional derivations of the classical result. The 
first derivation is based on a physically appealing master equation approach, in the spirit 
of the Boltzmann equation. The second derivation is a straightforward application of the 
Kubo formula. The calculation is done for G 1 and can be easily modified in order to get G 2 . 
The advantage of the Kubo formula approach is that it can be generalized to the quantum 
mechanical case, and it allows the incorporation of non-adiabatic and environmental effects. 



4. Classical derivation using a master equation 

We consider a ring with two scatterers: a moving scatterer go whose velocity is X, and a fixed 
scatterer gi. A collision of a particle with the moving scatterer implies that its velocity is 
changed v h-> v±2X, where the sign depends on whether the collision is from the right or from 
the left. The associated change in the kinetic energy is E i-> E + 2mvX + 0(X 2 ) respectively. 
There are two regions (x < and x > 0) on the two sides of the go scatterer. Accordingly we 
have four distribution functions that satisfy the following balance equations: 

o — * 

-gf = -lp7 v ]+8o\P-v] + (l-go){p7v] E - 2mvX (16) 
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^ = -b-V]+£o[prv] + (l-So)b-V] (18) 

The zero order solution in X is to have all the four distribution functions equal to some 
arbitrary function f(E). In the presence of driving, assuming that the system has reached 
a steady state, we still have to satisfy the two X-free equations, leading to 

px = gip-+a-gi)p+ (20) 

P- = giP+ + (1 - gi)P- (21) 
Substitution into the two other equations leads after linearization to 

P -(E) - p-(E) = -2mvxf l ~% ) ^ (22) 
\ga + gi-2gogi) oE 

and for the current we get 



7 = I l^- pr)ev= I %^-P. 

Jo l^\go + 



)ev (23) 

o ^ 

(1 -go)g\ 



g\ - 2gog\ 



T« <24 » 



With the assumption of zero temperature Fermi occupation this gives the cited result for G 1 . 
5. Classical derivation using the Kubo formula 

The generalized fluctuation-dissipation version of the Kubo formula (see Ref. 1 4 1 and further 
references therein) relates the generalized conductance to the the cross correlation function 
of the current I and the generalized force T = -cfHjdX. If X is the displacement X\ of the 
scatterer then 

r = -—=X 2 6'(x-X 1 ) (25) 
oX\ 

For the sake of comparison with previous results we assume zero temperature Fermi 
occupation. Then the Kubo formula takes the form 

{I(r)T{0))dT = ^-{QF) (26) 

where Q{E) = L/(7rhv F ) is the density of states. This density of states is proportional to the 
total "volume" of the network which is L. In the second expression we got rid of the time by 
introducing the notation 

<2= I(T)dT (27) 



It should be clear that both the generalized force T and the transported charge Q are functions 
in phase space, and that (...) stands for phase space average over position and velocity. For T 
we already have an explicit expression Eq.j25ll. Now we have to figure out what is Q. 

On the ring there are two scatterers, and one point x = xq where the current is measured. 
Hence the ring is divided into 3 segments. In addition, there are two possible directions 
of motion (clockwise, anticlockwise). Hence the phase space is divided into 6 regions. It 
is obvious that the outcome from Eq.j27t depends merely on which region the classical 
trajectory had started its journey in. In fact we need to consider only the 4 regions where 
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the particle starts in the vicinity of the moving scatterer, else T vanishes. So we have the " + " 
region between the moving scatterer and xq, and the " - " region on the other side between the 
two scatterers. Accordingly the four possible outcomes from EaA21t are: 

1 



fir 

Q- 

Q- 



2(1 -g T ) 
1 



2(1 



-gT) 

go 



gi(l - go) 



l-(l-gi)(l-go) 
gi-go+ gogi 



i-a-goa-go) 



er = 



go ~g\ + gog\ 
go + g\ ~ g0g\ 



q: 



(28) 
(29) 
(30) 
(31) 



go + gi -gogi 

The derivation of the above expressions is as follows. It is simplest if the particle starts in 
the " + " region, because then we can regard the two scatterers as one effective scatterer 
gr- Assume that at time t = the particle approach x = xq from the left. The charge that 
goes through the section after a round trip is suppressed by a factor (2gr - 1) due to the 
scattering (we sum the clockwise and the anticlockwise contributions). Thus we find that the 
total charge that goes through the section due to multiple reflections is a geometric sum that 
leads to Ea. (l28> . If we start in the " + " region in the opposite direction, then we have the 
same sequence but with the opposite sign and without the first term. Hence we get Ea. (l29> . 
Next assume that at t = the particle starts in the " - " region, and approaches go from the 
left. Then we can have at a later time a positive pulse of current. The probability for that is 
the geometric summation over go((l - go)(l - gi)) in ' C8C1 . Otherwise, we get a negative pulse 
of current, with a complementary probability that can be regarded as a geometric summation 
over gi((l - go)(l - gi)) inlcg "(l - go)- Thus the total current through the section, taking into 
account all subsequent multiple reflections (rounds) is given by Ea. J30> . A similar calculation 
leads to Ea.dSTl. 

Since there are only four possible values for Q the calculation of the phase space average 
becomes trivial: 

















l rdr . 


^ + 2l 


I™. 











Q- 



The integral over T is taken either within the " + " or within the " - " region. It is trivially 
related to the momentum impact and yields the result 



I 



Tdr — +mv F 



(32) 



Putting everything together we get the desired result for G l . With some minor modifications 
we can calculate G 2 using the same procedure. 



6. The quantum mechanical picture 

The Kubo formula holds also in the quantum mechanical case. But now I and T are operators, 
so it is more convenient to express the Kubo formula using their matrix elements. After some 
algebra one obtains the result: 

2Mm[X, lm ]^ m „ 



(E m - E n ) 2 + (r/2) 2 



For more details see Ref . 1 4 1 and further references therein. In the above formula it is assumed 
that only one energy level («) is occupied. If we have zero temperature Fermi occupation, then 
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we have to sum over all the occupied levels. The Kubo formula incorporates a parameter F 
that reflects either the non-adiabaticity of the driving, or environmentally induced "memory 
loss" due to decoherence. For a strictly isolated system in the strict adiabatic limit we have 
F = 0. Then we identify G as an element of Berry's field B, as explained in the introduction. 
The effect of T on B will be discussed below. 

We would like to see how the classical result can emerge in some limit from the above 
quantum expression. It turns out that this does not require a detailed calculation. We can use 
some topological properties of B in order to figure out the answer! The main observations 
that we further explain below are: 

(1) The B is divergence-less with the exception of Dirac monopoles 

(2) The monopoles are arranged in X space as vertical chains 

(3) The far field of B is like a two-dimensional electrostatic problem 

(4) Only non-compensated chains give net contribution 

As long as the occupied level n does not have a degeneracy with a nearby level, B is finite 
and divergence-less. Only at degeneracies can it become singular. It can be argued that these 
singularities must have their charge quantized in units of H/2 else the Berry phase would be ill 
defined. We have defined X3 — O as the Aharonov-Bohm flux through the ring. This means 
that if we change X? by Inh/e then by gauge invariance we have another degeneracy. This 
means that the Dirac monopoles are arranged as vertical chains, and that the average charge 
per unit length is e/(4n). Thus the far field of a Dirac chain is as in a two dimensional 
electrostatic problem. If we calculate the line integral of Eq.Q then we get, within the 
framework of the far field approximation, Q = 1. Thus we conclude that if we have several 
Dirac chains of the same "sign", then Q simply counts how many are encircled. 

We have to notice that if we have Fermi occupation, then the net contribution comes 
only from degeneracies of the last occupied level with the first unoccupied level. This is what 
we meant above (item 4) by "non-compensated". In order to avoid misunderstanding of the 
"compensation" issue let us discuss with some more details what happens if two neighboring 
levels n and m are occupied. With the level n we associate a field B ( "\ while with m we 
associate a field B (m) . In general B (m) + -B^K If we are near a degeneracy than we may 
say that B (n) emerges from a Dirac chain which is associated with level n, while B (m) emerges 
from a Dirac chain which is associated with level m. By inspection of Eq.j33>. taking into 
account that Im[i~„ m ] = -lm[I mn ], we realize that the two Dirac chains have opposite charge. 
Their corresponding fields do not cancel each other, but the total field is no longer singular, 
implying that the net charge is zero. 

In the quantum stirring problem we shall see that the X\ distance between non- 
compensated chains is simply half the De-Broglie wavelength Ag = 2n/kE- From this it 
follows that the amount of charge which is pushed by a very "large" scatterer is 

Qae ^± =e h xAXi (34) 
A E /2 n 

What happens if the cycle is not in the "far field" but rather passes through the distribution 
of the monopoles? To be more specific let us consider what happens to Q if we displace the 
scatterer a distance AXi. What is the dependence on X% \ Do we get the classical result as in 
Fig. 5? Obviously, in order to get the classical result the distribution cr(X l ,X2) should be in 
accordance with Eq.(ll5>. Strictly speaking this is not the case because we have a discrete set 
of monopoles rather than a smooth distribution of "magnetic charge". Still we can hope that 
cr(Xi,X2) would be classical-like upon course graining. We discuss further this issue in the 
next paragraphs. 
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If we make a pumping cycle in the vicinity of a monopole then it is obvious that the result 
would be very different from the classical prediction. What we expect to get in the quantum 
mechanical case is illustrated in the upper panel of Fig. 7. For a cycle that goes very close 
to a monopole the charge can be huge. In reality it is very difficult to satisfy the adiabatic 
condition near a degeneracy, or else there are always environmental effects. Either way, once 
we have a finite F, the result that we get for Q is smoothed. 

If the pumping cycle passes through a distribution of many monopoles then what we 
expect to get (as we deform or shift the cycle) are huge fluctuations as illustrated in the lower 
panel of Fig. 7. Again, the effect of either non-adiabaticity or environmental effects is to 
smooth away these fluctuations. The interested reader can find some further discussion of this 
point including a numerical example in 1121 . 



near 





Fig.7. Several pumping cycles are indicated in the left panels: It is implicit that each segment is closed as in Fig. 3. 
The black points represent degeneracies. For each pumping cycle one can calculated Q. The qualitative expectation 
for the outcome is illustrated in the right panels. In the upper illustration we assume that the pumping cycle encircles 
only one degeneracy, while in the lower illustration we assume that it encircles N degeneracies. In a later section we 
display numerical results that support the illustrated expectations. 

Coming back to the quantum-classical correspondence (QCC) issue, we realize that 
at best QCC can be satisfied in a statistical sense. So we ask whether the coarse grained 
cr(Xi,X2) agrees with the classical expectation Ea.( ll5l >. The answer which we give in the 
following sections, is that QCC is not realized in the case of simple non-chaotic models. In 
the "simple" cases we get a non-classical cr(X\ , X 2 ) and hence a different dependence of Q on 
X 2 . 
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7. The degeneracies in X space 

We can use the scattering approach in order to find the energy levels of a ring. In this approach 
the ring is opened at some arbitrary point and the S matrix of the open segment is specified. 
It is more convenient to use the row-swapped matrix, such that the transmission amplitudes 
are along the diagonal: 



SiE:X l ,X 2) = c-[ : ^-_ in / Vl7J , 



The periodic boundary conditions imply the following secular equation 

det(S(£;X 1 ,X 2 )-l) = (36) 
Using 

det(S - /) = det(S) - trace(S) + 1 (37) 
det(S) = (e'>) 2 (38) 
trace(S) = 2V?e iy cos0 (39) 

we get 

cos(y(E)) = Vi(^)cos(0) (40) 

In order to find the eigen-energies we plot both sides as a function of E. The left hand side 
oscillates between -1 and +1, while the right hand side may have a smaller amplitude. It is 
not difficult to realize that the only way to have two eigen-energies coincide is to get 

{0 = mod(27r) ) ( (p = n mod(27r) ) 

g = 1 or I g = 1 (41) 

7 = n„ cn n J [ y = n M n J 

where n is either even or odd integer that can be exploited (if we keep track over y) as a level 
counter. 

Both g and y depend on (E; Xi,X2). Since we want g to be maximal the condition for 
having a degeneracy involves 4 rather than 3 equations as we are going to see below. An 
immediate conclusion is that we have two types of Dirac chains: those that have monopoles 
in the plane of the pumping cycle (X3 = O = 0), and the others that have monopoles off the 
plane of the pumping cycle. 

In our model system we have two scatterers. One is the moving scatterer and the other is 
the rest of the network. The two are connected by arms of length La - X\ and Lb + X\ . The 
constants L A and L B can be absorbed into the definition of the surrounding network. Each 
scatterer is fully characterized by the set of parameters {gi, y it or,-, 0,}. Note that we do not 
absorb X\ into the definition of ao. After some algebra we find the following expressions for 
the transmission coefficient and for the total phase shift: 

8 - g ° gl (42) 

2 - go - gi + gogi + 2 y(l - g Q )(l - gi)cos(y + J\ + a + <*i - 2k E Xi) 

7 = 70+ 7i (43) 

where kg is the wavenumber that corresponds to the energy E. Thus the conditions for having 
a degeneracy take the form 



X3 = integer flux 

go(Xi) = gi 

a + a \ - 2k E X\ - n mod(27r) 
70 + 71 = n cvc „7r 



Xt, - half integer flux 

go{Xi) = gi m 
ao + a 1 - 2k E Xi - mod(27r) 

70 + 71 = "uddTT 
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We have highlighted the dependence on the parameters {X\,Xi,Xs). There is of course also 
an implicit dependence of {g h y,, a,} on the energy E. The conditions that are listed above are 
very intuitive: The system should have time reversal symmetry; The barriers should "balance" 
each other; The phases which are associated with the reflections should lead to destructive 
interference; And the total phase shift should respect the periodic boundary conditions. 

From Eq. J44b ~) we see that in general the X\ distance between degeneracies that belong 
to the same level is roughly half the De-Broglie wavelength as stated previously. The question 
that we would like to address is how these degeneracies are distributed with respect to X2- 



8. Quantum stirring in simple rings 

We would like to find the distribution of degeneracies with respect to Xi in the simplest model: 
a ring with two delta scatterers (see Fig.l). The arms that connect the two scatterers are of 
length La + X\ and Lb~X\. For the S matrix that represents the fixed scatterer (including the 
arms) we have 

-l 

(45) 

yi(E) = k E (L A + L B ) - arctan|^-X 2 J (46) 

ai (E) = k E (L A -L B ) (47) 

Since the dependence of go and g[ on the barrier "size" has the same functional form, the 
condition Eg. (144b) implies Xi — V irrespective of E. Thus we get that all the degeneracies 
are concentrated at the same Xi- This is clearly very different from the classically expected 
distribution. 

In Fig. 8 we display an example. The degeneracies that are associated with the first 7 
levels are indicated. Filled circles stand for = degeneracies, while hollow circles stand for 
(p = n degeneracies. Only the last (7th) level contributes non-compensated monopoles. The 
X\ distance between the non-compensated monopoles is roughly half De-Broglie wavelength. 

In Fig. 9 we show what happens to the degeneracies if we add a second fixed scatterer. 
We have chosen an additional scatterer that can be treated as a perturbation. The calculation 
was done using perturbation theory. We shall not present the details of this lengthy 
calculation here. For larger perturbations (not presented) we had to solve the secular equation 
numerically. This was done using an efficient algorithm 1 19 1. In any case, the purpose of Fig. 9 
is merely to demonstrate that once the symmetry of the system is broken the degeneracies 
spread out in the Xi direction. 
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Fig.8. The degeneracies in the double delta model of Fig. 1. We set La = 10.23 and Lg = 0, so that X[ measures 
the distance from the fixed scatterer. The "size" of the fixed delta scatterer is V = 1274.56. We use units such that 
m = h = 1. We assume that only the lower 7 levels are occupied. The filled circles are degeneracies on the flux 
zero plane and the empty circles are degeneracies on the flux n plane. The left graph shows the actual arrangement 
in the (Xi,X2) plane. Namely, all the degeneracies are on the line X2 = V. In the right graph the degeneracies were 
displaced for the sake of clarity. Only the 7th occupied level contributes non-compensated monopoles. 
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Fig.9. The degeneracies in the triple delta model of Fig. 1. Namely, to the model of Fig.8 we have added a delta barrier 
of "size" Vp = 10~ 5 , located at x = 7.61. This additional delta barrier can be treated as a small perturbation. As a 
result of this perturbation the degeneracies shift and spread out in the X2 direction. Degeneracies that belong to the 
same level are connected by a line. As in the previous figure only the 7th occupied level contributes non-compensated 
monopoles. 

The distribution (r(Xi,X2) in the case of a ring with a single fixed scatterer is very 
different from the classical prediction. Consequently also Q comes out very different from 
Eq.(ll5> [and see also Fig. 6]. The reader might be curious to know how Q depends on 
the "size" {X{) of the scatterer in the case of Fermi occupation. So we have calculated G 
numerically using Ea.J33>. and integrated over it to get Q. The numerical results are displayed 
in Fig. 10. Further analysis of the crossover from "near field" to "far field" cycles will be 
published in a separate work 1171 . 
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Fig.10. Several pumping routes are displayed in the left panel. For each of them Q has been calculated numerically. 
The results are displayed in the right panel. Note the agreement with the qualitative expectation that has been 
expressed in Fig. 7. The calculation is done for the double delta model of Fig. 1 with La = 1000.23 and Lg = 0. The 
"size" of the fixed barrier is V = 810.56. The energy level involved are n = 998 and m = 999. We use units such that 
m = h= 1. 



9. Quantum stirring in chaotic rings 

We would like to find the distribution of degeneracies with respect to Xi in case of a chaotic 
network (see an example in Fig. 1). Let us try to extend the approach that has been used in the 
previous section. A hypothetical illustration of gi(E) in the chaotic case is displayed in Fig. 1 1 . 
The universal conductance fluctuations of gi are characterized by a one parameter probability 
distribution P(g\;gi) which we discuss below. This probability distribution depends on one 
parameter, which we choose to be the average transmission^!. 




°0 200 400 600 800 

E 



Fig.ll. A hypothetical illustration of g[ (E) in the case of a complex "chaotic" barrier. Such a barrier can be modeled 
as a network (Fig. la), or it can be characterized using random matrix theory. The smooth curves are the transmission 
go(E;X2) of the delta scatterer for 3 different values of Xi- 

In order to get a degeneracy, a necessary but insufficient condition is that the transmission 
of the two barriers is equal (go(E;X2) = gi(E)). The solution of this equation can be 
determined graphically via Fig.ll. In fact in most practical applications we can assume that 
our interest is restricted to some small energy window such that the smooth E dependence 
of go can be neglected. So the equation is in fact #0(^2) = g\(E). For a given E we can find 
an such that this equation is satisfied. By playing with X\ we can satisfy the a related 
phase condition for having a degeneracy. But we still have to satisfy also the y related phase 
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condition, which leads to the quantization of the energy E. Hence the erratic XT* is sampled. 
Still it is reasonable to assume that the the distribution of the so-obtained X 2 values is not 
affected by this random-like sampling. We therefore conclude the following relation: 



Prob 



X 2 < Xf <X 2 + dX 2 ] = Prob[g (X2) < gi < go(X 2 + dX 2 )\ 



(48) 



This implies a simple relation between cr(X\ , X 2 ) and the probability function P{g\ ; g\) 

(r{X x , X 2 ) = const x ^®P( g0 (X2)) (49) 

UA 2 

Thus the problem of finding cr(X l , X 2 ) has reduced to the problem of finding P{g\ ; g\). 

We can now proceed in three directions: (A) To determine PQ from simple heuristic 
quantum chaos considerations; (B) To determine PQ from formal random matrix theory 
considerations; (C) To use reverse engineering in order to determine what is PQ that would 
give the classical result. It should be clear that universality can be expected only if £1 <K 1. In 
Fig. 12 we make a comparison between the outcomes of these three procedures for £1 = 0.001 . 
In the following paragraph we give the details of the calculation. 




0.06 



Fig. 12. A plot of the distribution P(gi ; gi) according to several different expressions. In this calculation we assume 
that the average transmission is gi = 0.001, which is represented in the figure by a vertical dashed line. The 
"heuristic" result is based on sampling of the random variable gi = gir][T]2 where r\ is Porter-Thomas distributed. 
The "RMT" result is based on Eg. 1501 . The "classical" result is based on Eg. 1511 . 

The heuristic approach is based on the idea that the transmission via a chaotic 
network depends on the amplitudes of the wavefunctions at the entrance and exit points. 
One might expect g\ = g\T]\7] 2 , where 77 has the Porter-Thomas distribution |20| 
PaoEin) = (1/ yllnrfjor^ 1 . This leads to the "heuristic" result in Fig. 12. In fact this result 
should not be taken too seriously. The formal RMT calculation 1211 of the probability 
distribution P(gv, gi) leads to the following expressions: 



PRMj(gi;gi) 



(2ln'g x )g\ 112 forg! «(£,) 2 « 1 
,(4£i/7r 2 )g- 3/2 for(g 1 ) 2 «g 1 « 1 

The small g\ approximation is universal: it merely assumes that the system has time reversal 
symmetry. It has been confirmed [22 1 that this universal behavior holds also for network 
systems. But for larger values of g\ there are deviations that has to do with semiclassical 
considerations. It is therefore in the latter region where one might expect quantum-classical 
correspondence. 



(50) 
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The probability distribution P(gugi) that would reproduce the classical result Eg.(fT3l 
via Eq.j49> is: 

Pa.(gi;gi)= , , ./'^'. c m (51) 



(gi+g?-2g lg p 

,1 ~ n Mn, 



with gy m 0.12gi. In order to compare with the RMT result we note that 



Pc(gUgl)~{ gfg -2 f 0I g?« gl «l (52) 

We see that in the large g\ region, where one might expect quantum-classical correspondence, 
there is no agreement between P C l() and P RMY Q. We suspect that P RM t() cannot be trusted there. 
Otherwise we have to conclude that Eq.j49> fails to take into account strong correlations in 
the arrangement of Dirac monopoles. Either way it seems that RMT alone is not enough in 
order to reproduce the classical result. 



10. The emergence of the classical limit 

With simple minded RMT reasoning we have failed to get a quantitative correspondence with 
the classical result. We therefore look for a different way to get an estimate for either B2 or 
cr(Xi,X2) in the case of a chaotic network. One obvious way is to use the result of Ref. [23 1 
regarding the distribution of degeneracies (diabolic points). The perturbation term which is 
associated with X2 is 

W=^=S(x-X 1 ) (53) 

OA 2 

and the density of the degeneracies should be [23 1 

(t{X u X 2 ) = ^g(E) 2 RMS[T„ m ] RMS[0V„ m ] <* RMS [TV™] (54) 

where Q(E) is the density of states. In the first equality it is implicit that the root mean square 
(RMS) of near diagonal matrix elements should be estimated. In fact only RMS[1V„ m ] is 
required in order to find the X2 dependence. For a quantum chaos system with time reversal 
symmetry the variance of the near diagonal elements equals half the variance of the diagonal 
elements |24|, leading to the second expression. 

There is a well known semiclassical recipe I25I I26I for calculating the variance of 
the near diagonal matrix elements r W nm . One should find the classical correlation function 
C(r) = <nV(f)nV(0)> - <0V> 2 , and then integrate over t. If "W were the current operator then 
{ r W) would be equal to zero, and we could proceed as in section 5. But in case of Ea. (l53> 
there is a problem: The sign of 'W(f) does not fluctuate, and it is essential to take into account 
the distribution of the delay times inside the network. Therefore there is no obvious relation 
to the transmissions g and g\. 

An optional possibility is to try to evaluate RMS['W / „„], where = \if/ bimi J 2 is the 
"intensity" of the wavefunction at the location of the scatterer. Obviously the result depends 
on both go and g\, and requires considerations which are at least as difficult as estimating 
universal conductance fluctuations. So it seems that we would run into the same problems as 
in the previous section. 

Still there is the option to calculate G l = B2 from the Green function of the system. 
This has been done in 1 12|: Writing the Green function as a sum over trajectories, we have 
expressed G l as a double sum over paths. If this double sum is averaged over the energy one 
obtains the diagonal approximation, leading to the classical result. At first glance the energy 
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averaging is not quite legitimate, because the energy is quantized. But one can justify this 
procedure in the case of a "quantum chaos system". We have further supported this claim by 
the numerical analysis of the chaotic network of Fig. 1 1 12|. We therefore conclude that for a 
chaotic network the distribution of degeneracies should be in accordance with Eq.( I15> . 

11. Conclusions 

As we translate a scatterer of "size" X 2 a distance AX\ along a single mode wire, the amount 
of charge which is pushed is 

Q = r(X 2 ) x —k F x AX] (55) 
n 

where k F is the Fermi momentum. If the scatterer is very "large" (X2 — > °°) then we expect 
to have r(X 2 ) = 1. This expectation is based on the "snow plow" picture that has been 
explained in the conclusion of section 2. This result is also confirmed by the formal BPT 
based calculation in the case of an open geometry. It also can be formally derived for a closed 
geometry using the "Dirac chains picture". In the latter case the key observation is that the X\ 
distance between contributing degeneracies is roughly half the De-Broglie wavelength. See 
Eq.lG!- 

Next we ask what happens to r(X 2 ) as X 2 becomes smaller. In the case of an open 
geometry the intuitive naive guess, which is based on the "snow plow" picture, turns out to be 
correct. Namely, r(X 2 ) — 1 - go is simply the reflection coefficient: Some of particles are not 
"pushed" by the scatterer because of its partial transparency. In the case of a closed geometry 
we have shown that the classical result for r(X 2 ) is modified: now it depends also on the 
overall transmission of the device. See Ea. (ll3> . 

It is important to realize that the classical result for r(X 2 ) is in complete agreement with 
the common sense expectation. Namely, we have < r(X 2 ) < 1, and the dependence on the 
"size" of the scatterer is monotonic. But once we go to the quantum mechanical analysis we 
have a surprise. The results that we get are counter-intuitive. They are most puzzling (Fig. 10) 
in the case of the simplest model, in which the ring contains only one fixed delta barrier (V). 
As we decrease X 2 the transported charge Q becomes larger(!). Moreover, once X 2 becomes 
smaller than V, the coefficient r(X 2 ) changes sign. This means that as we push the particles 
"forward" the current is induced "backwards". 

The reason for the failure of our intuition is our tendency to regard "adiabatic transport" 
as a zero order adiabatic approximation, while in fact it is based on a first order analysis (for 
a detailed discussion see section 4 of |17|). As a parameter in the system is changed, the 
induced current can be in either direction. 

In order to understand the route towards quantum-classical correspondence it is essential 
to figure out how the degeneracies spread out in X space. As the system becomes more 
complex, we get for r(X 2 ) a result that resembles the classically implied one. The resemblance 
is at best only on a coarse grained scale: the quantum result has strong fluctuations. These are 
related to universal conductance fluctuations. 

We have made an attempt to deduce from RMT considerations the "chaotic" distribution 
of the degeneracies, and hence the dependence of r(X 2 ) on X 2 . The quantitative results do 
not agree. We therefore suspect that RMT considerations alone are not enough in order to 
establish quantum-classical correspondence. Rather we had used 1121 semiclassical tools in 
order to establish this correspondence. 
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